On the Kinematic Origin of the Luminosity-Pulse Lag 
Relationship in Gamma-ray Bursts 
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ABSTRACT 

This paper presents an interpretation based on gamma-ray burst source kinematics for the 
relationship found by Norris et al. between peak luminosity and energy-dependent pulse lag. 
I argue that the correlation should instead be between number luminosity and pulse lag. This 
interpretation improves the least-squares fit of this correlation for the known bursts by 25 percent 
or more. It also suggests a distance estimation scheme. I propose that this relationship is due to 
the variation in line-of-sight velocity among bursts. This interpretation allows one to speculate 
on the range of gamma-ray burst expansion velocities or the size of their jet opening angles. 

Subject headings: gamma rays: bursts — gamma rays: theory 



1. Introduction 

It has recently been reported by Norris et al. 
(2000a, b) that there is a relationship between 
the peak luminosity of gamma-ray bursts (GRB) 
and the pulse time lag between BATSE energy 
channels. Specifically, they find that the peak 
luminosity L p k is anticorrelated with the dura- 
tion of the lag, At, according to L p k ~ 1.3 x 
(At/0.01 sec) -115 xlO 53 ergs/sec . This intriguing 
relationship is useful because it probes the prop- 
erties of an ensemble of bursts and thereby can 
provide clues of the global dynamics of GRBs. It 
also has the potential for distance estimation to 
GRBs. 

Several workers have attempted to explain the 
origin of the pulse lag, or so called pulse paradigm 
(Norris et al. 1996). Daigne & Mochkovitch (1998) 
were able to produce a pulse lag by modeling the 
burst as a series of internal shocks, but obtained 
timescales an order of magnitude too large. How- 
ever, as Wu & Fenimore (2000) point out, the 
lag timescale is much longer than the synchrotron 
cooling time. Panaitescu & Meszaros (1998) mod- 
eled the kinematics of internal shocks and found 
no evidence of pulse peak lag. As of yet, the source 
of the pulse paradigm remains a mystery. This pa- 
per does not attempt to explain the origin of the 



pulse lag, but assumes that it derives from a pro- 
cess common to all bursts and has some proper 
decay timescale At' in the reference frame of the 
emitter. 

In this paper I will present an interpretation 
of the luminosity-lag relationship as being due to 
relativistic motion of the emitting region toward 
the observer. Thus it is a purely kinematic effect. 

2. Photon Number Luminosity and Cos- 
mology 

In Norris et al. (2000a) a peak 7-ray luminosity 
Lpk was defined from the observed peak number 
flux F p k, a luminosity distance Dl derived from 
the observed redshift z and a given cosmology, and 
a mean emitted photon energy e, constant for all 
bursts: 



L pk = 4ireF pk D 2 L 



(1) 

However, one does not know the intrinsic energy 
of the photons emanating from the pulses. There- 
fore multiplication of all GRB fluxes by an av- 
erage photon energy e, while useful to estimate 
GRB peak luminosites, obfuscates the underlying 
relationships in the dynamics of GRBs. In other 
words, photon number flux F p k is the observed 
quantity and thus photon number luminosity N p k 
can be the only calculable luminosity. 
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The luminosity distance Dl is denned so that 
redshift of photons is accounted for. The photon 
number luminosity does not depend on photon en- 
ergy redshift and thus is given by 



N, 



pk 



^F pk Dj 
(1 + *) 



(2) 



In Figure 1 is shown the peak number lumi- 
nosity N p k versus spectral lag At for six bursts 
with known rcdshifts. I fit the data by minimiz- 
ing Q 2 = S(log(7Vp fe . mo dci)-log(iVp fejda t a )) 2 where 
Npfc.modd = A(At)- p . The best fit for lags CCF31 
0.1 (see Fig. 1) is 



N pk = 8.6 x 10 56 At" 



0.98 



photons sec 



(3) 



with Q 2 — 1.38. By comparison, a similar fit 
of L pk from Eqn (1) versus At reproduces the 
Norris et al. (2000a) result L pk oc At" 1 - 15 with 
Q 2 = 1.73. For CCF31 0.5, N pk oc At" 89 with 
Q 2 = 0.45 contrasts with L pk oc At -103 with 
Q 2 = 0.79. With six datapoints, these are small 
number statistics to be sure. However, the im- 
provement in correlation N pk vs. At over that of 
L pk vs. At via reductions of Q 2 by 25% or greater 
is clear. 

3. Relativistic Kinematics 

The relation N pk oc At" 98 , where 0.98 is al- 
most 1, is a strikingly simple statistical property 
for an ensemble of GRBs to obey. In fact this rela- 
tionship suggests a simple explanation. I propose 
that this proportionality is a result of variations 
in line-of-sight relativistic velocities of GRBs. 

For a burst expanding with Lorentz factor 7 = 
l/\/l — (i 2 , (3 = v/c and at redshift z, the ob- 
served number luminosity N bs (photons sec -1 ) 
varies as 



N obs = 



7 7V' 



(4) 



where TV' is the proper isotropic number luminos- 
ity for a non-expanding source. This expression 
is valid for emission that appears isotropic to the 
observer (i.e. if the emission is in the form of a 
jet, then the opening angle 8 > I/7, where I/7 
is the relativistic beaming angle). Also, as men- 
tioned previously, if one assumes that the spectral 
pulse lag is due to some proper decay timescale 
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Fig. 1. — Peak photon number luminosity N pk 
versus spectral pulse lag for six bursts with known 
redshifts. Spectral cross-correlation function lag 
between BATSE channels 3 and 1 (CCF31) for 
regions down to 0.5 and 0.1 of peak intensity were 
obtained from Norris et al. (2000a). The line of 
best fit for 0.1 (squares) is oc At -0 98 . 



At', possibly due to cooling or deceleration, then 
we observe a lag 



At o6s = 7(1 - /3cos0)(l + z) At' 



(1 + z)At' 



2 7 



(5) 



for 7 3> 1, where the angular dependence is re- 
moved by the strong effect of attenuation of num- 
ber flux F received from emission moving at in- 
clination angles greater than 1/7; F ~ (7(1 — 
/?cos0)(l + z))~ 3 - Thus one can focus one's at- 
tention only on emitters moving directly toward 
the observer. Combining these then yields the in- 
verse relationship 



N obs oc At 



obs 



(6) 



This relationship is redshift independent and thus 
is valid for luminosity N pk = (l + z)N O0S and spec- 
tral lag At = At b s I {1 + z) as would be observed 
local to the burst. Thus we have an explanation of 
the relation shown in Figure 1. I propose that the 
variety in the spectral lags and number luminosi- 
ties of bursts depend solely upon the relativistic 
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motion of the emitter along the line of sight. Cer- 
tainly one expects variation among bursts to cause 
deviations from this relationship, however the cur- 
rent data suggests this effect to be dominant. 

There are two possibilities for the nature of this 
variation (see Section 6). First, we are observing 
an intrinsic range of burst expansion speeds, rang- 
ing from fast bursts having little or no spectral 
lag to slower bursts having longer lags. The sec- 
ond possibility is that emitting material is in the 
form of a jet with some angular dependence on the 
emitter's Lorentz factor 7(6*) where the Lorentz 
factor would have a peak at the centroid of the jet 
j(8 = 0) = -f m ax and monotonically decrease on 
some angular scale greater than 1 /j for increasing 
inclination angles 8. In either case we observe only 
the emission dominated by those emitters moving 
with trajectories inclined within an angle 1/7 from 
our line-of-sight. 

4. Redshift Estimation 

The above discussion allows the estimation of 
redshift z from the observed peak number flux F pk 
and the observed spectral lag Ar j s . Fenimore & 
Ramirez-Ruiz (2000) also estimate redshifts using 
GRB variability to indicate luminosity. From Eqn 
(2) one has 



given observed peak number flux and spectral lag 



N pk =4irF pk (l + zfD 2 A 



(7) 



where Da(z) is the angular diameter distance (e.g. 
Hogg 1999) which is a function of the cosmology 
and redshift z, but does not vary much over the 
interval 1 < z < 3. For example, here I follow Nor- 
ris et al. (2000a) and choose (Q m ,^a) = (0.3,0.7), 
for which Da ~ OADh to within 10% over this 
range of z, where Dh = c/Hq is the Hubble Dis- 
tance (choose H = 65 km sec -1 Mpc -1 ). Also, an 
observed spectral lag At bs varies from the local 
spectral lag At by 



At = 



At obs 
(1 + *) 



(8) 



One can put these components into the N pk (At) 
curve of bursts with known redshifts 



pk 



A 

At ' 



(9) 



where A = 8.6 x 10 56 photons as taken from the fit 
in Figure 1. Thus one can estimate a redshift for a 



z ~ 



4nF pk D 2 A At obs 



- 1 ~ 



1.5 



y/F pk: At bs 



- 1 



(10) 



where flux F pk is measured in photons cm -2 sec -1 
and At bs is in seconds. It is worth noting that a 
well determined cosmology gives a well determined 
Da (z) and thus a refined redshift z. Conversely, 
many independently determined redshifts might 
give cosmological information. However, the pri- 
mary source of error in this relation is not the un- 
certainty in Da{z), but in the deviations of burst 
luminosities and lags from the relation of Eqn (9) 
as can be seen in Figure 1. 

5. Observed Peak Pulse Energy 

In order for there to be a lag At of the arrival of 
a pulse over an energy range fixed with respect to 
the observer (e.g. the range between BATSE chan- 
nels 1 and 3) that varies inversely with Lorentz 
factor (Eqn 5), the peak energy of the pulse must 
decay exponentially as 



E pk ~ E' pk o jexp 



J7 



(ii) 



where E' pk is the peak energy in the emitter 
frame, t is observer time, 7 is the Lorentz factor of 
the emitter and T pk is the proper decay timescale 
in the emitter frame. Thus we infer an exponential 
decay law for the pulse peak energy. 

It is important to note that this exponential 
evolution can be very rapid (< 10 msec) (Norris 
et al. 1996, 2000a) and is distinct from the evo- 
lution of the pulse envelope (the so-called FRED) 
that has been studied by many authors (e.g. Rydc 
& Svensson 2000; Liang & Kargatis 1996; Norris 
et al. 1996; Fenimore & Bloom 1995). 

One outstanding question in the study of GRB 
pulses is the apparent invariance among bursts of 
the peak energy E mx of the E 2 N E flux spectrum 
(E mx ~ a few x 100 keV) where is number flux 
per photon energy. One might expect this value to 
vary with burst luminosity, but it remains roughly 
constant (to a factor of a few) over a wide variety 
of (long) GRBs. 

If the pulse photon number flux Ng (which, for 
demonstration, we assume to be peaked at E pk ) 
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increases as the peak photon energy E p k of Eqn 
(11) decreases, then there will be a maximum in 
E 2 pk NE- This may happen as follows: the nascent 
pulse emitter is extremely hot and compact. As it 
evolves, it cools (E p k drops) and it expands or be- 
comes progressively more optically thin, releasing 
an ever increasing photon rate Ne- This is anal- 
ogous to early-time supernova lightcurves (East- 
man et al. 1994). In this picture, uncovering the 
flux dependence Nsit, 7) will give insight into why 
E mx ~ constant. For instance Ne oc i 1 " 7 = j lnt 
with Eqn (11) demonstrates this behavior. 

6. Discussion 

In this letter I have argued that the correlation 
of Norris et al. (2000a) is between number lumi- 
nosity and spectral lag and I interpret this as being 
due to the variety, among bursts, of relativistic ve- 
locities at which emitting regions move toward the 
observer. Suppositions as to the physical emission 
process have been deliberately avoided in order to 
highlight this purely kinematic effect. 

There are two possible scenarios by which this 
might happen. The first is simply that the varia- 
tion is due to the variation of expansion Lorcntz 
factor among bursts. From Eqn (5), this implies 
that the ratio of maximum to minimum spectral 
lags is the same as the ratio of maximum to minu- 
mum Lorentz factors 

Imax At max (-IO\ 

■y ■ ~ At ■ ~ ^ > 

where observed lags are seen to range over roughly 
two orders of magnitude. Thus if the fastest bursts 
have j m ax ~ 100, then the minimum bursts would 
have jmin ~ 1 and thus may be only mildly rela- 
tivistic. 

The second possibility is that GRB ejecta is di- 
rected in a jet such that the Lorentz factor ^(9) 
has some maximum -f(9 = 0) = jmax and mono- 
tonically decreases with increasing angle. Thus 
the function 7(0) will determine the relative num- 
bers of bursts observed with given spectral lags. 
Band (1997) has noted that the distribution of 
bursts is strongly peaked at small lags indicat- 
ing that the Lorentz factor of emitting ejecta, 
7(0 < 6o/2) ~ j m ax, is basically constant over 
the jet opening angle 9 , with a fairly narrow edge 
region 9 > 9 of decreasing 7. Thus let us estimate 



a characteristic opening angle 9q by assuming that 
7 = Imax is constant for 9 < 9 /2 and decreases to 
order unity within a narrow boundary or "edge" 
region < 9 — 9^/2 < l/^ max . Thus the solid 
angle subtended by "face" bursts (i.e. observed 
at inclinations 9 < 9 /2) is 7r(# /2) 2 , and that for 
"edge" bursts is ir9o/ , ~fmax and the ratio is: 

"edge" bursts _ 4 
"face" bursts -y ma xOo 

From Fig. 3 of Norris et al. (2000a), one can es- 
timate this ratio, by defining short lags < 0.1 sec 
as "face" and the longer lags as "edge" bursts, to 
be ~ 2/5. Assuming j m ax = 100, then this gives 
5° < 6> < 10°. With more statistics of the distri- 
bution of bursts along the curve of Eqn. (6) one 
can fit the function 7(0) much more precisely. 

This work was performed under the auspices of 
the U.S. Department of Energy by University of 
California Lawrence Livermore National Labora- 
tory under contract W-7405-ENG-48. 
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